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Abstract 

An analysis is given of the chromoelectric dipole moment of quarks and of the neutron 
in an MSSM extension where the matter sector contains an extra vectorlike generation of 
quarks and mirror quarks. The analysis includes contributions to the CEDM from the ex¬ 
change of the W and the Z bosons, from the exchange of charginos and neutralinos and the 
gluino. Their contribution to the EDM of quarks is investigated. The interference between 
the MSSM sector and the new sector with vectorlike quarks is investigated. It is shown that 
inclusion of the vectorlike quarks can modify the quark EDMs in a significant way. Further, 
this interference also provides a probe of the vectorlike quark sector. These results are of 
interest as in the future measurements on the neutron EDM could see an improvement up 
to two orders of magnitude over the current experimental limits and provide an instrument 
for a further probe of new physics beyond the standard model. 
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1 Introduction 


New sources of CP violation beyond those that exist in the Standard Model are needed 
to explain baryogenesis and are also worthy of study in their own right as possible probes 
of beyond the standard model physics (for reviews see e.g., [H [21 [3l 0]). Such sources can 
also induce electric dipole moment in elementary particles which can be signihcantly larger 
than those expected in the standard model Bi. In this work we are specihcally interested 
in the electric dipole moment (EDM) of the quarks arising from the chromoelectric dipole 
operator. Thus the electroweak sector of the standard model produces an EDM which is 
10“^° ecm mum and it lies beyond the possibility of its observation in the foreseeable 
future. As mentioned in particle physics models beyond the standard model it is possible 
to generate much larger values for the EDM. In this work we focus on one such model - an 
extension of the minimal supersymmetric standard model (MSSM) with a vectorlike multi- 
plet [S]. Such an extension is anomaly free and thus the nice quantum properties of MSSM 
are maintained. Further, vectorlike multiplets arise in a variety of settings such as in grand 
unihed models and in string and D brane models [El 001 HI]. Vectorlike generations have 
been considered by several authors since their discovery would constitute new physics (see, 
e.g., 0[ini[II102l[ISllia05l[ISlEllISl[ISlEni[2Ill22]). Such models have new sources of CP 
violation and thus can generate substantial size dipole operators. For that reason they are 
interesting models to consider in the context of lepton and quark EDMs. In a recent work 
we analyzed the electric dipole operator in such a setting [23] and in this work we analyze 
the chromoelectric dipole operator in the extended MSSM model and its contribution to the 
electric dipole moments. 

Before discussing the EDM in the new class of models, it is relevant to recall the situation 
regarding the lepton and quark EDMs in MSSM. Here it is well known that MSSM has a 
SUSY CP problem, i.e., that the EDM predicted with SUSY phases 0(1) are typically in 
excess of the experimental upper limits. A number of remedies have been offered in the 
past to remedy this problem. These include a hne tuning of the phases to be small [23], 
suppression of the EDM by large sparticle masses [25], suppression of the EDM where various 
contributions conspire to cancel, i.e., the cancellation mechanism [261 EZ] as well as other 
possible remedies (see, e.g., [2E])- It has also been suggested that the EDM be used as 

a probe of new physics beyond the standard model [221 |30l jSTl [181 [22|- Specihcally the 
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experimental limits on the EDMs can be used as vehicles to probe a new physics regime not 
accessible otherwise to current and future detectors. 

The outline of the rest of the paper is as follows: In section 2 we give a brief description 
of the model and describe the nature of mixing between the vector generation and the 
standard three generations of quarks. In section 3.1 we discuss the loop contributions to 
the chromoelectric dipole moment of the up quark and the down quark that arise from the 
exchange of the W boson in the loop. In section 3.2 we give an analysis similar to that of 
section 3.1 for the exchange of the Z boson. In section 3.3 we compute the contribution from 
the exchange of charginos in the loop and in section 3.4 a similar analysis for the exchange of 
neutralinos in the loop is given. Finally in section 3.5 we give the analysis for the exchange 
of the gluino in the loop. In section 4 we discuss the method for the computation of the 
neutron dipole moment using the quark dipole moments. In section 5 we give a detailed 
numerical analysis of the contributions to the quark CEDM and to the neutron CEDM for 
a variety of parameter points in the extended MSSM model. Here we also discuss the use of 
the neutron EDM as a probe of high mass scales. Conclusions are given in section 6. Further 
details of the calculational aspects of the analysis are given in sections (7-9). 

2 The Model 

Here we briefly describe the model and further details are given in the appendix. The 
model we consider is an extension of MSSM with an additional vectorlike multiplet. Like 
MSSM the vectorlike extension is free of anomalies and as discussed in section 1 vectorlike 
multiplets appear in a variety of settings which include grand unihed models, string and D 
brane models. Here we focus on the quark sector where the vectorlike multiplet consists of a 
fourth generation of quarks and their mirror quarks. Thus the quark sector of the extended 
MSSM model is given by Eq. (1) and Eq. (2) where, 

(tx) ^ * = 1,2,3,4. (1) 

0=.(«|)~(3-.2.-t) ; T,~(3 .i4) ; (2) 

The numbers in the braces show the properties under SU{3)c x SU{2)l x U{1)y where the 
hrst two entries label the representations for SU{3)c and SU{2)l and the last one gives 
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the value of the hypercharge normalized so that Q = + Y. We allow the mixing of 

the vectorlike generation with the hrst three generations. Specihcally we will focus on the 
mixings of the mirrors in the vectorlike generation with the first three generations. Details 
of these mixings are given in Eq. (43). Here we display some relevant features. In the up 
quark sector we choose a basis as follows 


^R ^R ^R ^4/?) , — (fi Tl Cl Ul tn) ■ (3) 

and we write the mass term so that 

- £“ = S(M„)a + hx., (4) 

The interaction of Eq. (43) lead to the up-quark mass matrix which is given by 
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This mass matrix is 

not hermitian and a 

bi-unitary transformation is needed to diagonalize 


it. Thus one has 

D“'^(M„)T)^ = diag(m„,,m„2,m„3,m„4,m„5). (6) 


Under the bi-unitary transformations the basis vectors transform so that 
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A similar analysis can be carried out for the down quarks. Here we choose the basis set 

Vr = (pR Br Sr cIr hiR^ , rjp = (pL Br Sr cIr b^R) . (8) 

In this basis the down quark mass terms are given by 

- C‘P = fR{Md)iqL + \i.c., (9) 
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where using the interactions of Eq. (43), has the following form 


(yxvxl\f2 
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In general /i 3 , /i 4 , /is, / 13 , / 14 , / 15 , / 13 , / 14 , / 15 , /le, /^r, /^s can be complex and we define their phases 
so that 

hk = |/lfc|e*^^ /I'fc = |/i'fc|e*^S h'l = \h'l\e^^'^ . (11) 

The sqnark sector of the model contains a variety of terms including F -type, D-type and 
SUSY soft breaking terms. The details of these contributions to sqnark mass square matrices 
are discussed in section [71 

3 The analysis of Chromoelectric Dipole Moment Op¬ 
erator 

The chromoelectric dipole moment dP is the coefficient of the effective dimension 5 operator 
which is defined by 

.C; = (12) 

where is the gluon field strength and T“ are the SU{3) generators. The quarks will 
have five different contributions to the CEDM arising from the W, Z, gluino, chargino and 
neutralino exchanges. We denote these contributions by d^{W), d^{Z), d^{g), d^{x~^) and 
We discuss each of these contributions below. 

3.1 W exchange contribution to quark CEDM 

For the up quark the W- exchange contribution arises from the left diagram of Fig. (1) using 
the interaction of Eq. (13), i.e.. 


C^wu = Wl + G%^PR]d, + h.c.. 


i=l j=l 


(13) 
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Figure 1: W and Z exchange contributions to the CEDM of the up quark. Similar exchange 
contributions exist for the CEDM of the down quark where u and d are interchanged and 
W~^ is replaced by W~ in the diagrams above. 


where GY and GY are dehned in section 8. The contribution of the W-exchange graph to 
is given by 


= lib E 

167r^ ^ mw V ruw rriu 






where /i(ri,r 2 ) is a form factor given by 


^ , (4 + ri - r2)x - Ax'^ 

h{ri,r2)= dx—— --—^ 

Jo 1 + Yi - r2 - l)x + r2X^ 


In the limit when r 2 is very small as the case here, this integral gives the closed form 

T / r.\ 2 r 1 1 2 3ri Inri 1 

^'(’'‘.0) = [i + jn + -r, + J ■ (1 

The W contribution to the down quark CEDM is given by 

<'?(«-) = ^ t J) . (1 


3.2 Z exchange contribution to quark CEDM 

For the Z boson exchange the interactions that enter with the up type quarks are given by 


= Zp E E fipVKjiPL + 


j=l i=i 
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where the couplings and are dehned in section 8. Using this interaction the com¬ 
putation of the Z exchange contributions to the up quarks is given by the loop diagram to 
the right in Fig. (1). Its contribution is 


5 

9s 

i=l 

For the Z boson exchange, the interactions that enter with the down type quarks are 
given by 



^Im(U 

mi 


uZ ^uZ*\ T 

J-'l 


m: 
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m: 


W4 


mi 


(19) 


5 5 

-Cmz = Z,Y^ + CgPRlrf., (20) 

j=l i=l 

where the couplings and Cf^ are as dehned in section 8. A calculation similar to that 
of the up quark CDEM gives a contribution to the d-quark moment so that 
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9s 

1671-2 2 ^ 
i=l 

3.3 Chargino exchange contribution to CEDM 

In this section we discuss the interactions in the mass diagonal basis involving squarks, 
charginos and quarks. Thus we have 

5 2 10 

= E E Z + Cfi^PR)r'dt + hx., ( 22 ) 

7 = 1 i=l fc=l 

and 



5im(c£c£-)/i 

11L fy 


2 ’ 2 


( 21 ) 


5 2 10 

= E E E d^CfiPL + C’^,PR)rut + hx., (23) 

7=1 i=l k=l 

where the couplings (2^“, (2^“, and and are as dehned in section 8. The loop 
contributions to the up -quark CEDM arise from the right diagram of Fig. (2). Their 
contribution to CEDM of quarks using Eq. (22) and Eq. (23) are given by 
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Figure 2: Left diagram: Supersymmetric loop contributions to the CEDM of the up-quark 
from the diagram involving the exchange of neutralinos and up-squarks. Right diagram: 
Chargino and down-squark loop contribution to the CEDM of the up quark. Similar loop 
contributions exist for the CEDM of the down quark, where u and d are interchanged, u and 
d are interchanged and is replaced by in the diagrams above. 
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where h{ri,r 2 ) is given by 


'‘I 2 

X — X 


= / dx - 7 

Jq 1 (ri - r 2 - l)x + r 2 X^ 


(26) 


In the limit when r 2 is very small as is the case here we have the closed form 


^ , 1 / 2ri In ri 

4(D, 0) = -7^ ( 1 + ri + 


2 (ri - 1)2 


1 - ri 


(27) 


3.4 Neutralino exchange contribution to CEDM 

We now discuss the interactions in the mass diagonal basis involving up quarks, up squarks 
and neutralinos. Thus we have. 
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' ^u—u—yX 


+ <mPR)XiUk + hx., 


(28) 


i=l j=l k=l 
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The interaction of the down quarks, down squarks and neutralinos is given by 


4 4 10 

= E E E + CllkPR)x°A + li'C'. (29) 

i=l j=l k=l 

where the couplings C'^ and C'^ as given in section 8. Using the interactions of Eq. (28) 
the neutralino exchange contribution to the CEDM of the up-quark is given by 
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(30) 


Similarly using the interactions of Eq. (29) the CEDM of the down quark is given by 
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3.5 Gluino exchange contribution to CEDM 


(31) 


3 3 8 5 10 

-^qqg = ^ ^ ^ ^ ^ [Cli^Pl + CR^^PR]gaqti + h-c., (32) 

ji=l k=l a=l 1=1 m=l 

where the couplings 0^^ and Cri^ are dehned in section 8. Using Eq. (32) the gluino 
exchange contribution to the up quark CEDM arising from the loop diagrams of Fig. 3 is 
given by 
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Similarly using Eq. (32) the gluino contribution to the down quark CEDM is given by 


(33) 



9s^s 
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si 
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Here and are given by 




Figure 3: Left diagram: Supersymmetric loop contributions to the CEDM of the up-quark 
arising from the exchange of gluino and up squarks with the gluon emission from the internal 
up squark line. Right diagram: Same as left diagram except that the gluon emission is from 
the internal gluino line. Similar loop contributions exist for the CEDM of the down quark, 
where u and d are interchanged, u and d are interchanged. 


- D-nA - OSA “ KiAA", (35) 


and 


= (ofA + i^l'A + ^I'A + KiA - . 

where / 5 (ri,r 2 ) is the loop function defined by 

x + 8x‘^ 

h{ri,r2)= dx-— -——- 

Jq 1 (ri - r2 - l)x + r2X^ 

In the limit where r 2 is very small as is the case here we get the closed form 
T f ^ An on i 2rilnri ISlnri 

0) = 7777-777 1 lOri -2Q + 
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(36) 


(37) 


(38) 


4 The neutron CEDM 


As discussed in the previous section, the total contribution to CEDM of the quarks consists of 
hve contributions arising from the exchange of the W, the Z, the charginos, the neutralinos, 
and the gluino, so that 


^ = <(«')+<if(Z)+<if(C)+ <if(x“) + <f(9). 9 = «.<*■ (39) 
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The contribution of the chromoelectric operator to the EDMs of quarks can be computed 
using dimensional analysis [31] . The contribution to the quark EDM arising from is given 
by 

< = ( 40 ) 

where is approximately equal to 3.4. The factor t]*" brings the electric dipole moment from 
the electroweak scale down to the hadronic scale where it can be compared with experiment. 
To obtain the contribution to the neutron EDM from the quark EDM, we use the non- 
relativistic SU{6) quark model which gives 




(41) 


5 Numerical analysis of neutron EDM 

The current experimental limit on the EDM of the neutron is [35] 

\dn\ < 2.9 X 10"^® ecm (90% CL). (42) 

It is expected that a higher sensitivity by as much as two orders of magnitude more sensitive 
than the current limit may be achievable in the future [3B] . 


We present now a numerical analysis of the neutron CEDM Erst for the case of MSSM and 
next for the MSSM extension. The first analysis involves no mixing with the mirror genera¬ 
tion and with the fourth sequential generation and the only CP phases that appear are those 
from the MSSM sector. Thus in this case all the mixing parameters, given in Eq. (11), are 
set to zero. The second analysis is for the MSSM extension where the mixings of the mirror 
generation and of the fourth sequential generation with the three generations are switched 
on. The results are given in Table and Figs. |4 IT In the analysis, in the squark sector 
we assume mf = M? = M? = M? = M? and mf = M? = Ml = M? = M% = M? = 
= M^j. = M? . To simplify the numerical analysis further we assume mg = mg = mo. 


bs' 


At 


Additionally the trilinear couplings are chosen so that: = At = At = Ac = A^ = A 

and Aq = Ah = Ab = As = Ad = A^h- The input parameters are such that the sparticle 
spectrum that enters the loop are consistent with the current experimental limits from the 
LHC in each of the cases, i.e., with or without mixing. 
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Contribution 

( 

:i) 

(ii; 

) 

Up 

Down 

Up 

Down 

Chargino, d^ 

2.49 X 10-29 

-1.29 X 10-26 

2.16 X 10-29 

-2.08 X 10-26 

Nentralino, d^° 

-2.49 X 10-32 

4.75 X 10-29 

-2.90 X 10-32 

5.47 X 10-29 

Glnino, d^ 

3.42 X 10-29 

-4.24 X 10-23 

7.49 X 10-23 

2.06 X 10-26 

Total, dq 

5.90 X 10-29 

-1.32 X 10-26 

7.71 X 10-23 

-1.42 X 10-28 

EDM, d® 

-2.70 

X 10-26 

-6.83 X 10-28 

Chargino, d^'(x^) 

-3.41 X 10-3° 

-2.15 X 10-27 

-2.89 X 10-39 

-3.40 X 10-27 

Nentralino, d^(x°) 

-4.54 X 10-32 

-1.73 X 10-23 

-5.30 X 10-32 

-2.00 X 10-23 

Glnino, dg{g) 

5.51 X 10-29 

1.37 X 10-27 

1.21 X 10-27 

-6.63 X 10-26 

Total, dg 

1.40 X 10-29 

-2.58 X 10-23 

3.26 X 10-23 

-1.89 X 10-26 

CEDM, d^ 

-3.49 

X 10-23 

-2.53 X 

10-26 


Table 1: An exhibition of the chargino, nentralino and glnino exchange contribntions to the 
quark and the neutron EDM, CEDM and their sum for the case when there is no mixing 
of the vectorlike generation with the three generations. The analysis is for two benchmark 
points (i) and (ii). Benchmark (i): 6^ = 3.3 x 10“^, ,^3 = 1 x 10“^. Benchmark (ii); 
9fj_ = 4.7 X 10“^, .^3 = 3.6. The common parameters are: tan/? = 40, mo = mg = mg = 3000, 
|mi| = 185, Imal = 220, |A“| = 680, \A^\ = 600, |/i| = 400, rUg = 1000, |h 3 | = \h'^\ = \h'^\ = 

\h 4 \ = \K\ = \K\ = \h\ = |h'I = \h''\ = \h\ = M = IhsI = 0, 6 = 2 X 10-2, 6 = 2 X 10-3, 

aA^ = 2 X 10“^, aAd = 3. All masses are in GeV, all phases in rad and the electric dipole 
moment in ecm. 

We discuss now in further detail the cases without and with mixing with the vectorlike 
generation. We begin with the case with no mixing. In table we give the individual con¬ 
tributions to the up and down quark EDM and CEDM, namely, the chargino, the nentralino 
and the glnino contributions. The W and Z contributions are not shown since they are absent 
in this case of no mixing with the vectorlike generation and the fourth sequential generation. 
The different contributions are given for two benchmark points (i) and (ii), where in the 
hrst, the neutron EDM dominates the neutron CEDM and in the second, the opposite is the 
case. The chargino and glnino contributions are the main contributors, while the nentralino 
contribution is suppressed. Note that the total neutron EDM, \dn\, obtained by adding 
and d^ in the table satisfy Eq. (42). Another observation is the largeness of the down quark 
contribution in comparison with its up quark counterpart. This is attributed to the large 
value of tan /? which tends to enhance the down quark couplings. 
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( 

:i) 

(ii; 

) 

Contribution 

Up 

Down 

Up 

Down 

Chargino, 

7.65 X 10-30 

-6.91 X 10-27 

7.08 X 10-30 

-8.27 X 10-27 

Nentralino, 

3.93 X 10-32 

6.90 X 10-30 

3.91 X 10-32 

7.32 X 10-30 

Glnino, d^ 

-2.01 X 10-23 

-5.35 X 10-27 

2.25 X 10-28 

5.91 X 10-27 

W Boson, df 

3.77 X 10-30 

3.46 X 10-28 

3.77 X 10-30 

3.46 X 10-28 

Z Boson, dg 

8.02 X 10-31 

3.05 X 10-20 

8.02 X 10-33 

3.05 X 10-20 

Total, dq 

-1.88 X 10-28 

-1.19 X 10-20 

2.37 X 10-28 

-1.97 X 10-27 

EDM, df 

-2.41 

X 10-20 

-4.13 X 10-27 

Chargino, d^'(x^) 

-7.66 X 10-31 

-8.98 X 10-28 

-6.95 X 10-33 

-1.07 X 10-27 

Nentralino, d^(x°) 

7.17 X 10-32 

-2.52 X 10-20 

7.14 X 10-32 

-2.68 X 10-20 

Glnino, dg{g) 

-4.57 X 10-28 

2.44 X 10-20 

5.13 X 10-28 

-2.70 X 10-20 

W Boson, d^{W) 

-2.97 X 10-30 

2.29 X 10-28 

-2.97 X 10-30 

2.29 X 10-28 

Z Boson, dg{Z) 

1.46 X 10-30 

-1.11 X 10-28 

1.46 X 10-30 

-1.11 X 10-28 

Total, dg 

-1.24 X 10-28 

6.38 X 10-27 

1.38 X 10-28 

-7.56 X 10-27 

CEDM, d^ 

8.54 X 10-27 

-1.01 X 10-26 


Table 2: An exhibition of the chargino, nentralino, glnino, W and Z exchange contribntions 
to the qnark and the nentron EDM, CEDM and their sum for the case when there is mixing 
of the vectorlike generation with the three generations. The analysis is for two benchmark 
points (i) and (ii). Benchmark (i): 6 ^^ = 4x 10“^, ^3 = 1.12. Benchmark (ii): 9^ = 4.6x 10“^, 
,^3 = 4.71. The common parameters are: tan/? = 40, mo = m^ = ttiq = 5500, \mi\ = 185, 
|m2| = 220, |Ag| = 680, |Aq| = 600, |/i| = 400, rUg = 1100, m^ = 300, m^ = 240, = 320, 

m 4 b = 280, |h 3 | = 1.58, \h'^\ = 6.34 x lO-^, |h"| = 1.97 x lO-^, |h 4 | = 4.42, \h'^\ = 5.07, 
|h"| = 12.87, Ihl = 6.6, Ih'sl = 2.67, \h'^\ = 1.86xl0-^ {hel = 1000, \h7\ = 1000, Ihgl = 1000, 
6 = 2 X 10-2, 6 = 2 X 10-3, = 2 X 10-2, = 3 , X3 = 2 X IO- 2 , y' = 1 x 10-^, 

X'' = 4X 10-3, X 4 = 7X 10-3, x' = x" = 1X 10-3, X 5 = 9 x 10-3, x '5 = 5 X 10-3, x' 5 ' = 2 x 10-3, 
Xe = X 7 = Xs = 5 X 10-3. All masses are in GeV, all phases in rad and the electric dipole 
moment in ecm. 

Next we consider the case with mixings. The results are presented in table for two 
benchmark points (i) and (ii). Here in addition to the chargino, the nentralino and the 
glnino exchanges one also has W and Z exchanges. The analysis shows the dominance of the 
EDM over the CEDM for benchmark (i) while opposite is the case for benchmark (ii). The 
total EDM for each benchmark point satisfies the experimental constraints of Eq. (42). Here 
we note that the EDM and CEDM are constrained not only by the experimental limits on the 
MSSM spectrum, but also by the limits on new quarks. Thus for the benchmarks presented 
in Table 2, the vectorlike quarks have masses gotten by diagonalization of the matrices of 
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Eq. (5) and Eq. (10) and are given in TableThe results of Tableare consistent with [57] . 
More stringent constraints on these masses will be available at the LHC RUN-II. 


Mirror Up Quark mf = 1037 GeV 

Mirror Down Quark mb' = 740 GeV 
Fourth G Up Quark = 1057 GeV 

Fourth G Down Quark = 1260 GeV 


Table 3: An exhibition of the masses of the heavy extra quarks corresponding to the param¬ 
eter space of table 



Figure 4: Variation of neutron GEDM \d^\ versus Mx {Mx = Ih^l = Ih-^l = Ihgl), for four 
values of mo- From top to bottom at Mx = 4 TeV, mo = mg = m^ = 2.0, 2.5, 3.0, 3.5 TeV. 
Other parameters have the values tan/3 = 14, |mi| = 185, |m 2 | = 220, |/i| = 350, |Ao| = 680, 
|Aol = 600, mT = 300, ms = 260, mg = 1000, m^ = 320, m^b = 280, lhal = 1.58, 
Ih'l = |h"| = RMx, \h,\ = 4.42, |h'I = \h'l\ = RMx, \h\ = 6.6, |h'| = \h''\ = RMx, 
R=lx 10-3, Og = 3.98, 6 = 6 = 4.52, 6 = 2.42, = 5.0, = 1.14, xs = 2.38, 

X' = 4.92, X 3 = 2.58, X4 = 4.86, x '4 = 1-6, xl = 1-37, Xs = 1-14, X 5 = 4-39, Xs = 2.38, 
Xe = 4.92, X 7 = 2.58, xs = 4.86. All masses are in GeV and phases in rad. 


Next we give an analysis of the quark GEDMs dependence on the mass scales as well 
as on the GP phases both in the MSSM sector as well as the new sector. Thus the GEDM 
depends on the mass scale of the vectorlike sector and in the MSSM sector it depends on 
the universal scalar mass mg, and on the gaugino mass scales. Further, it has dependence 
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on several CP phases both from the MSSM sector as well as from the vectorlike sector. We 
discuss the dependence of the CEDM on the mass scales hrst, and specihcally on the mass 
scales Mx (from the vectorlike sector) and on mo and on irtg. 

Fig. 1^ gives the dependence of the effect of the vectorlike generation on CEDM where 
we exhibit CEDM vs Mx, where Mx = \h^\ = \h^\ = Ihg] and that \h'^\ = \h'^\ = \h^\ = 
Wil = 1 ^ 5 ! = 1 ^ 5 ^ = RMx while R = 1 x 10 “^. We note that the allowed range of values 
for R is highly constrained. Thus smaller values of R will not produce interesting results 
while larger values are likely to upset the quark masses for the hrst three generations. The 
analysis shows that CEDM lower than the upper limit can be obtained and masses in the 
TeV range may be probed using the constraint given by Eq. (42) which should undergo 
further rehnements in the future. The curve corresponding to mo = 2 TeV is characterized 
by a dip at Mx ~ 1.9 TeV. This dip quickly widens and is replaced by a shallow drop for 
mo = 2.5 TeV and then disappears completely for larger values of rriQ. The variation of the 
CEDM eventually levels off for higher values of Mx and mo. Further analysis shows that 
the dip is caused by a sudden drop in the mass of the lightest up squark mass for Mx ~ 1.9 
TeV in this region of the parameter space. The analysis of the dip is rather involved but 
arises as a result of the competition among the different components of the chromoelectric 
dipole moment operators, i.e, the W, Z, chargino, neutralino and gluino contributions. The 
analysis of Fig. makes clear the very sensitive dependence of the CEDM on the vectorlike 
mass scale and exploration of this dependence is one of the primary motivations of this work. 
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Figure 5: Variation of neutron CEDM \d^\ versus the scalar mass mg (mo = mg = mg), 
for five values of Mx, [Mx = Ihg] = \hj\ = Ihgl). From top to bottom at mo = 5 TeV, 
Mx = 1.5, 2.0, 2.5,3.0, 5 TeV. Other parameters have the values tan/3 = 14, |mi| = 185, 
|m2| = 220, |/i| = 350, \Aq\ = 680, |Aq| = 600, mr = 300, uib = 260, rUg = 1000, mu = 320, 
mu = 280, Ihsl = 1.58, \h'^\ = Ih'^] = RMx, \h4\ = 4.42, \h'^\ = \h'l\ = RMx, Ih^] = 6.6, 
|h'I = \h''\ = RMx, R=lx 10-3, Og = 3.8, = 6 = 4.52, = 2.42, = 5.0, = 1.14, 

X 3 = 2.38, X 3 = 4.92, Xs = 2.58, X 4 = 4.86, = 1-6, xl = 1-37, Xs = 1-14, Xs = 4.39, 

Xs = 2.38, Xe = 4.92, xi = 2.58, Xs = 4.86. All masses are in GeV and phases in rad. 


Another way for looking at Fig. |^is to plot the CEDM against mo for several values of 
Mx while R is fixed at 1 x lO-^ in the same region of parameter space. This is done in 
Fig. The plot shows peaks between 2 and 3 TeV and then the CEDM decreases gradu¬ 
ally for increasing values of mo. The peak is more pronounced for small values of Mx and 
disappears for larger values (for Mx = 5 TeV, here). The peaks occur in regions where mo 
and Mx are comparable in size as shown in this region. All values of the CEDM obtained 
in this region of the parameter space lie below the current upper limit. 

In Fig. I^we investigated the dependence of CEDM on the vectorlike mass Mx and found 
that there is a very significant dependence of the CEDM on Mx- It is of interest also to ex¬ 
amine if the EDM shows a similar dependence on Mx- In Fig. 6 we exhibit this dependence 
where \d^\ is plotted against Mx for the same set of mo values as in Fig. Again as in the 
case of CEDM we find that EDM also has a sensitive dependence on Mx- We note here that 
the analysis of this work for EDM is more general than the analysis of [23]. Thus in [23| 
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Figure 6: Variation of neutron EDM \d^\ versus Mx {Mx = \hQ\ = \hy\ = j/isl); for four 
values of tuq. From top to bottom at Mx = 1 TeV, mo = m^ = = 2.0, 2.5, 3.0, 3.5 TeV. 

Other parameters have the values tan/3 = 15, \mi\ = 185, |m 2 | = 220, |/i| = 350, |y4g| = 680, 

|y4g| = 600, ruT = 300, ms = 260, mg = 1000, m^ = 320, m^b = 280, [hsl = 1.58, 
j/i'gl = |/i"| = RMx, \h^\ = 4.42, \h'^\ = \h'i\ = RMx, \h\ = 6 . 6 , \h'^\ = \h'^\ = RMx, 

R = lx 10-3, Of, = 5 X 10-3, 6 = 2 X 10-2, 6 = 2 X 10-3, ^3 = 4.0, oa- = 2 x 10-^, 

= 3.0, X 3 = 2 X 10-2, Xs = 1 X 10-3, X 3 =^x 10-^ X 4 = 7 X 10-3, x '4 = x'i = 1 X 10-3, 
X 5 = 9 X 10-3, X 5 = 5 X 10-3, xl = 2 X 10-3, xe = X 7 = Xs = 5 x 10-3. All masses are in 
GeV and phases in rad. 


we considered only the mixings of the three generations with the mirror generation so that 
the quark matrices were 4x4 and the squark square matrices were 8x8 and the parameter 
Mx did not appear in that work. On the other hand, in this work we are considering mixing 
of the three generations with the full vectorlike generation consisting of the mirror and the 
sequential fourth generation. As a consequence the quark mass matrices are 5x5 and the 
squark mass squared matrices are 10 x 10 and this time we have also the dependence on the 
vectorlike mass Mx- Thus the analysis of this work is more general than of the work of |23] . 
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Figure 7: Variation of neutron CEDM versus the gluino mass, for four values of 
tan/3. From bottom to top at rrig = 5 TeV, tan/3 = 10, 20, 30, 40. Other parameters have 
the values \mi\ = 170, |m 2 | = 220, |/i| = 450, |y4g| = 680, [Agl = 600, uIq = rUg = 3700, 
mr = 300, tub = 260, mu = 320, = 280, = 1.58, jh'gl = 6.34 x lO'^, |h"| = 

1.97 X 10-2, \hi\ = 4.42, \h'^\ = 5.07, \hl\ = 2.87, Ih] = 6 . 6 , \h'^\ = 2.67, lh' 5 'l = 1.86 x 10-^ 
\h\ = \h 7 \ = Ihsl = 1000, Og = 2.6 x lO’^, 6 = 2 x lO'^, ^ = 2 x lO’^, ^3 = 1.6, 
ttAjf = 2 X 10-2, a^d = 3.0, X 3 = 2 X IO- 2 , ^3 = 1 x 10-^ Xs = 4 x 10-^ X 4 = 7 x 10-^ 
X 4 = X 4 = 1 X 10■^ X 5 = 9 X 10-3, X 5 = 5 X 10-3, xs = 2 X 10-3, xe = X? = Xs = 5 x 10-3. 
All masses are in GeV and phases in rad. 
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Figure 8: Variation of up and down quark contributions to the neutron CEDM for tan/? = 
40. Other parameters are the same as in Fig. 


Next we study the dependence of the CEDM on the gluino mass. This is given in Figs. 
and[^ Thus in Fig. the variation of the neutron CEDM, is plotted against the gluino 
mass, rUg. It is shown that CEDM values lower than the current experimental upper limit 
can be obtained in the given parameter space. The neutron CEDM decreases for increasing 
values of rUg, but eventually levels off at around zero for some values of tan/?. However, for 
other values of tan/?, (e.g. tan/? = 40), the neutron CEDM levels off but turns negative. 
This phenomenon can be understood by analyzing different contributions to the CEDM 
as shown in Fig. Specifically one finds that the negative contribution to the CEDM 
arises from the chargino exchange loop contribution, Since we are not applying any 

GUT constraints, the masses of the chargino and the gluino can be treated as independent 
parameters and thus as we increase the gluino mass, the chargino contribution remains 
unchanged and eventually dominates as the gluino mass gets large and makes the CEDM 
negative for rUg > 20 TeV. We note here in passing that the W and Z contributions in this 
region of the parameter space are negligible compared to the other exchange contributions. 
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Figure 9: Variation of neutron CEDM \d^\ versus for four values of tan/3. From bottom 
to top aX 9^ = 1 rad, tan/3 = 10, 20, 30, 40. Other parameters have the values \mi\ = 170, 
|m2| = 220, |/i| = 400, |Aq| = 680, |Aq| = 600, rug = tuq = 8000, rUg = 1000, rriT = 300, 
ms = 260, m 4 i = 320, m^b = 280, Ihg] = 1.58, Ih'gl = 6.34 x lO'^, \h'^\ = 1.97 x lO'^, 
|h 4 | = 4.42, \h'^\ = 5.07, \h'i\ = 2.87, j/isl = 6.6, \h'^\ = 2.67, |h"| = 1.86 x 10-^ Ihe] = Ihy] = 
\hs\ = 1000, 6 = 2 X 10-2, 6 = 2 X 10-3, es = 2.6, = 2 x 10-^, = 3.0, X 3 = 2 x 10-^, 

x'g = 1X 10-3, x's' = 4X 10-3, X4 = 7X 10-3, x', = x'l = lx 10-^ X5 = 9X 10-3, x's = 5X 10-3, 
Xs = 2 X 10-3, xe = X 7 = Xs = 5 X 10-3. All masses are in GeV and phases in rad. 


As discussed already, it is of interest to study the dependence of CEDM on the CP phases 
in the MSSM sector as well as in the new sector. Fig. shows the variation of the neutron 
CEDM versus 9^, the phase of /r. The CP phases are the source of the CEDM and the sensi¬ 
tivity that the CEDM shows in response to the variation of 9^ is obvious. The parameter /i 
appears in the chargino and the neutralino mass matrices. It exists also in the squark mass 
squared matrices, so one can see that the chargino, the neutralino and the gluino contribu¬ 
tions are affected by this parameter and its phase. The electroweak contributions, i.e, W and 
Z components are independent of the magnitude and the phase of /i. Values of \d^\ below 
the current upper limit can be obtained for several values of tan /3, whereas values above the 
limit appear for larger tan/3. 
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Figure 10: Variation of neutron CEDM \d^\ versus for four values of |/i 8 |. From bottom 
to top at Xe = 1 rad, j/igl = 1100, 1200, 1300, 1400 GeV. Other parameters have the values 
tan/3 = 34, |mi| = 185, \m 2 \ = 220, |/x| = 350, |2l^| = 680, = 600, = mi = 3600, 

mr = 300, ms = 260, mg = 4000, m 4 t = 320, m 4 b = 280, [hsl = 1.58, \h'^\ = 6.34 x 10“^, 

|h"| = 1.97 X 10-2, |h 4 | = 4.42, \h'^\ = 5.07, \h{\ = 2.87, Ih] = 6 . 6 , \h'^\ = 2.67, \hi\ = 

1.86 X 10-^ \he\ = Ihrl = 1100, 9^ = 0.1, = 2 x lO'^, ^2 = 2 x lO'^, ^3 = 3.6, 

ttA- = 2 X 10-2, a^d = 3.0, X 3 = 2 X IO- 2 , Xg = 1 x 10-^ Xs = 4 x 10-^ X 4 = 7 x 10-^ 

Xi = X 4 = 1 X 10-3, X 5 = 9 X 10-3, x's = 5 X 10-3, xi = 2 x 10-3, X 7 = X 8 = 5 X 10-3. All 
masses are in GeV and phases in rad. 


Next we investigate the dependence of GEDM on xe which explores a new sector of the 
theory as it is the GP phase that arises in interactions involving the mirror quarks and the 
fourth generation quarks. An analysis of the dependence of GEDM on xe is exhibited in 


Fig. 10 Aside from h^, other mass parameters that arise because of the new sector are hi 


and hg- The dependence of the GEDM on IhsI is also exhibited in Fig. 10 Quite remarkably 
GEDM is sensitive to both the mass scale and the phase that enters in the new sector. 
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Figure 11: Variation of the neutron EDM, \d^\ (solid curve), the neutron CEDM, \d^\ 
(dashed curve), and the total neutron EDM, (dotted curve), versus .^ 3 , the phase of 

the gluino mass, for tan/3 = 40. Other parameters have the values \mi\ = 185, \m 2 \ = 220, 
|/i| = 400, I^qI = 680, \Aq\ = 600, tuq = rug = 5000, tut = 300, tub = 260, rUg = 1500, 
rriit = 320, rriib = 280, jhal = 1.58, \h'^\ = 6.34 x lO'^, \h''\ = 1.97 x lO'^, \h^\ = 4.42, 
\h'^\ = 5.07, \h'l\ = 2.87, IhsI = 6 . 6 , \h'^\ = 2.67, lh' 5 'l = 1.86 x 10-^ Ihe] = Ihj] = j/isl = 1000, 
9^ = 4.7 X 10-^ 6 = 2 X 10-2, 6 = 2 X 10-^ ua- = 2 x lO'^, = 3.0, Xs = 2 x lO'^, 

x '3 = 1X 10-3, x' 3 ' = 4 X 10-3, X4 = 7 X 10-3, = X 4 = 1 X 10-3, X5 = 9 X 10-3, x's = 5 X 10-3, 

X 5 = 2 X 10-3, xe = Xt = Xs = 5 X 10-3. All masses are in GeV and phases in rad. 


Finally, it is of interest to look at the total electric dipole moment obtained by adding the 


electric and the chromoelectric dipole moments. Fig. 11 shows the variation of the EDM, the 


CEDM and the total EDM against the gluino phase, ^ 3 . The analysis of Fig. pT] shows that 
while the EDM may dominate the CEDM for some values of ^3 the opposite may happen for 
a different range of ^ 3 . The analysis also suggests constructive interference between EDM 
and CEDM in some parts of the parameter space while there is destructive interferences in 
other parts (i.e., for 0 < ^3 < tt) leading to the cancellations mechanism |261E7|. 


6 Conclusion 

In this work we have given an analysis of the chromoelectric dipole moment of quarks and 
of the neutron arising in an extension of MSSM where there is an additional vectorlike 
generation of quarks in the matter sector. Such an extension brings in new sources of CP 
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violation which can contribute to the chromoelectric dipole moment of quarks. The work 
presented here consists of analytical results on hve different types of contributions to the 
chromoelectric dipole moments of quarks which include both non-supersymmetric as well as 
supersymmetric loop contributions. In the non-supersymmetric sector we have contributions 
arising from the exchanges of the W and Z bosons in the loops, while in the supersymmetric 
sector we have exchanges involving charginos, neutralinos and the gluino in the loop. We 
have also carried out a detailed numerical analysis of their relative contributions. Specihcally 
it is found that there exists strong interference effects between the MSSM sector and the 
vectorlike quark sector which can drastically change both the sign and the magnitude of 
quark EDMs. We have also investigated the possibility that the neutron EDM can be used 
as probe of the TeV scale physics. These results are of import as future experiment can 
improve the current limits up to two orders of magnitude and thus the quark EDMs provide 
an important window to new physics beyond the standard model. 

Acknowledgments: PN’s research is supported in part by the NSF grant PHY-1314774. 


7 Appendix A: Squark mass matrices 

In this Appendix we give further details of the model discussed in section 2. As discussed 
in section 2 we allow for mixing between the vector generation and specihcally the mirrors 
and the standard three generations of quarks. We also allow for mixing between the mirror 
generation and the fourth sequential generation assuming R parity conservation (for a recent 
review of R parity see [33] )• The superpotential allowing such mixings is given by 


^1l] 


W = + y2H{Q‘^^fL + 

+ ysH\q^2L^2L + yz^2y2L^2L + + yi^^y^L^sL + y5Hiq{]^bh + y^Hiq\it'i 

+ h^eijQ^''q{]^ + q^^^ + + h'J)2iBL 

+ h'^il^TL + Khl^BL + Kil^TL + + hibl^BL + hdlLTi - , (43) 


Here the couplings are in general complex. Thus, for example, y is the complex Higgs mixing 
parameter so that y = |p|e®®'". The mass terms for the ups, mirror ups, downs and mirror 


downs arise from the term 


C 


1 d‘^W 

2 dAidAj 


h.c.. 


(44) 
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where ip and A stand for generic two-component fermion and scalar fields. After spontaneons 
breaking of the electroweak symmetry, {{H\) = Vi/\/2 and (hf|) = n 2 /v^), we have the 
following set of mass terms written in the fonr-component spinor notation so that 

= ^(M„)^L + hS(^rf)hL + h.c., (45) 


where the basis vectors are defined in Eq. (3) and Eq. (8). 

Next we consider the mixing of the down sqnarks and the charged mirror sdowns. The 
mass sqnared matrix of the sdown - mirror sdown comes from three sonrces: the F term, the 
D term of the potential and the soft SUSY breaking terms. Using the snperpotential of the 
mass terms arising from it after the breaking of the electroweak symmetry are given by the 
Lagrangian 

+ Aoft , (46) 

where Cp is dednced from Fi = dW/dAi, and —Cp = Vp = FiF* while the Cp is given by 
-Cd = cos^ dw cos 2(3{iLi*L - bpbl + cpcl - spsl + upup - dpdl + uJIl - hpblp 
+ ^rB*r - TrT^} + sin^ 9w cos2f3{-^iLil + ^IrFj^ - ^cpcl + ^crc}^ 

- \uluI + + ^TrT*^ - ^fpfl - ^ipbl - ‘^bRb*ji 

“ - ^dRdji + -BrBj^ 

+ ^BlBI - -UpilL + -^hnilR - ~ ^iRKR]- (47) 

For £soft we assume the following form 


-£.oft = + Mlqllqlp + 

+ MlbAjilp + Mli-LtlL + Mlbtpbl, + MliTpil, 

+ + Mlbtpbl^ + MlBlBp + Mlflfp 

+ eij{yiAbHlq{J)lL - y'^AtH^^dUiL + Vb^AbHldljilL - y'^AuH^^diJlL + yz^sHldiib^L 

- + y,A,H{qljbl^ - y',A^H^,qiJl^ + y^ApBlQ^^^TL - y'^ApH^.Q^^Bp + h.c.} 

(48) 


Here Mip, Mf, etc are the soft masses and At, Ab, etc are the trilinear couplings. The trilinear 
couplings are complex and we dehne their phases so that 


Ab=\Ab\F‘^^^ , At=\At\F^^^ ,••• . 


(49) 
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From these terms we construct the scalar mass squared matrices. Thus we dehne the scalar 
mass squared matrix M| in the basis {Sl, BL,bfi, Bji, sl, sr, dL, dR,h 4 L,b 4 ji). We label the 
matrix elements of these as {Mj)ij = M?- which is a hermitian matrix. We can diagonalize 
this hermitian mass squared matrix by the unitary transformation 

. (50) 

Similarly we write the mass squared matrix in the up squark sector in the basis 
iR,TR,CL,CR,UL,UR,i 4 L,i 4 R). Thus here we denote the up squark mass squared matrix in 
the form {M?)ij = mfj which is also a hermitian matrix. We can diagonalize this mass square 
matrix by the unitary transformation 

£>“'Af?£>“ = . (51) 

8 Appendix B: W, Z, x^, x^, g couplings with quarks 

8.1 W-quark -quark couplings 

The couplings that enter in the W -quark -squark interactions of Eq. (13) are dehned so that 

0%. = + DZ,Dtii + DtziDiu + (52) 

G"'. = (53) 

8.2 Z-quark -quark couplings 

The couplings that enter in the Z -up-quark interactions of Eq. (18) are dehned so that 


Cll = + DZiDui + DZ,Dlu + + V.DTyDl^l (54) 


and 


C'Z = + D^RtPRi. + DZjDlu + DZiD]m) + i.BS.Grk]. (55) 


COS 9 


w 


where 


xi = - - - sm 9w, yi = -- sm 9w 


(56) 
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The couplings that enter in the Z-down-quark interactions of Eq. (20) are dehned so that 


Cf_ = 


cos 9w 


(57) 


and 


Cf = 

JrC 'I 


[l/2(-DR5j-DR5i + + DR3jD%3i) + X2D'^2j^R2i\^ (^8) 


■\d* 


)d* 


d* 


■\d* 


cos 6 


w 


where 


-2. 

X2 = -^ + 2 Sin 9 


w, ^2 — 2 sin 9w- 


(59) 


8.3 Chargino-quark-squark couplings 

The couplings that enter in the chargino-up-quark-down-squark interactions of Eq. (23) are 
given by 


cjs; =»(-«:„ - i.,v-DZiDi, 

- + VlDHybit), (60) 
CS =a(->^XaDZ,bit - >‘.UaD'S,,bi - K,U,,Dl\,bi 

- KttUaDTiibiok - >‘tU,,DT,M 

+ + U„Dll^Di + U,,Dl\^bi, + U,iDllbl,), (61) 

The couplings that enter in the chargino-down-quark-up-squark interactions of Eq. (22) 
are given by 


cfi =3{->‘X;2 D%,d:;, - 


eft =3(-^uV,^D-fptt - KtV„_DttMk - ’^XaKjbi.t - 

- '^BV.^Dtyb'it + v,,Dtt,b:;t + Vt,D%bik + VnDtMk + KiKib^t), 

where 

{mT,mb,ms,md,m4b) 


{^HT 1 ^3si Ibid: ^4b) 


y/2mw COS /3 


(62) 

(63) 

(64) 


(^B, Kt: ^C1 ^U1 «4t) 


{mB:mt:mc:mu:mu) 
'/2m w sin / 


(65) 
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and 

= diag(m--, m--). (66) 

8.4 Neutralino-quark-squark couplings 

We first give the discuss the couplings that enter the the interactions in the mass diagonal 
basis involving up quarks, up squarks and neutralinos of Eq. (28). Here we have, 



26 



Similarly for the couplings that enter the the interactions in the mass diagonal basis 
involving down quarks, down squarks and neutralinos of Eq. (29) we have 

(75) 






d* T^d 


)d* f^d 


and 


elk - SkjDtMk + ^jOtuDik - S.iDtuDik + UbiDtMk - dkjDt.b 


d 

3k 


+ limDtbMk - hk,DtkkDt„k + aB,Dt2 A - 'lB,Dt2.Dtk), 


■\d* f^d 


d* f^d 


■\d* 


(76) 


where 


asj 

iBj 


gmsXlj _ 

2mw sin [5 ’ 

1 1 ^ sin^ Ow 

3 3 cos Ow 


and 


I^Bj 

Sbj 


^ w' , 3 

^ 

gmBXjj 

2mw sin /3 


w 


A-t (i 


t siii'^ 9„,- 


(77) 

(78) 


Oi4i)j 


gm^ijX^j 

2mw cos (3 ’ 


abj 


_ gm^bX*^ _ 

04bj „ n t ^bj 

2mw COS [5 

and where 


gmbX^j _ _ 

2mw cos 13 ’ 

grribXl- _ ^ ^ 

2mwcos/3 ’ 


grUsX: 


3i 


2mw cos [3 ’ 
grrisXl- _ 
2mw cos f3 ’ 


Oldj 


3dj — 


gmdXsj 

2mw cos (3 
gmdX*- 
2mw cos [3 


(79) 

(80) 


Pibj Pbj /3sj I3dj 

jAbj = Ibj = Isj = Idj = 

Here X' are dehned by 


.^eX'* H_ — 

3 cos 6 


-X. 


w 


2j 


+ ^sin^ 


1 V' , l^sin^^vv^/ 
3 “'^^ 3 cos»„, 


(81) 

(82) 


= Xu cos 9w + X 2 i sin 9w 
X2i = -Xu sin 9w + X 2 i cos 9w, 


(83) 

(84) 
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where X diagonalizes the neutralino mass matrix and is dehned by 


X'^M^oX = diag ym^o,m^o,m^o,m^oj . (85) 

8.5 G luino-quark-squar k-couplings 

The couplings that enter in the gluino-quark-squark interactions of Eq. (32) are given by 

= {dzA - - dzA - DZiA - dzAA'^. m 

and 

Ch,„ = {DTAn + DZA + DZA + DtnA - DZAA^ (87) 

where .^3 is the phase of the gluino mass. 


9 Appendix C: Mass squared matrices for the scalars 

We dehne the scalar mass squared matrix Mj in the basis (6^,5^, bn, Br, sr, sr, cIl, dR, b^i, & 4 r) 
We label the matrix elements of these as = Mf- where the elements of the matrix are 

given by 


MI 2 = M| + 


Ml^ = M| + 


2 


vl\y2\ 

2 

2 


vflyi 

2 

2 


V2\y2\ 

2 


+ \hi\ + \h'^\ + 1^11 + + o^z cos2/3sin 6 'ty, 


1 

3 ’ 


+ IhsI + \h'^\ + Ihgl + \he\ + ruz cos2/3 -sin Ow 


1 1 



^55 = ^II + + l^sT - COS2/3 Q - ^ sin^ 


Mm = + |/iir - 2^1 cos2/3sin^ Ow, 

Mji = + |/i3p - m| cos 2/5 Q - ^ sin^ ^h/ 

5^88 = Ml^ + + |/ilT - cos 2/3 sin^ • 

5^99 = + ^^2^^ + l^eT - COS2/3 Q - ^ sin^ 


'99 - ''''4L 


5^1010 - + 


"^lll/sp I I/, |2 5-2 0/0 • 2/1 

—- ^\h7\ --m2Cos2/3sm e'vi/ 

Zi o 
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= MTi = 


2* ^^21/2 ^3 , 2 


^ M^3 = Ml* = ^{v,Al - /it;2), Ml = Mil = 0, 

Ml = Mil = hlhl, Ml = Mil = 0, Ml = Mil = h'lhl, Ml = Ml = 0, Ml = MH = hlh, 
Mlo = Ml, = 0, Ml = Ml = 0, Ml = Ml = ^{v2A*^ - Ml = Ml = 


r2 _ n^2* _ n _ n^2* _ ^'2^31/2* , Viy^h'l ^ 


Ml = Ml = Q,Ml = Ml = 


l2 


+ 


, Ml — Ml — 0, 


72 


2 n^ 2 * '>Jihly5 , V 2 y 2 h(i ^_^2 .. 2 * 


Ml = Ml = 


+ 


l\yf2 _ 1^2* _ Q 


U3 


72 


+ 


72 


72 72 

; M^^ = M 53 = 0 , M^q = Mg 3 = , 


Ml = Ml = I '^ll/l^3 /\^2 _ T\/r2* _ n t\/t 2 _ T\/r2* _ u u'* 


Ml = Ml = 0, Ml = Ml = h.hl, 

M 39 = Mgg = 0 , M^,q = MgQ 3 = h4hj, 

Ml = Ml = 0, Ml = Ml = 


72 


72 ’ 


r2 _ A4^2* _ n A/r2 _ a j2* _ ^^27^7 , Vih'^yl 


Ml = Ml = 0, Ml = Ml = 


72 


+ 


72 ’ 

^2 _ J^2* _ 0 - M^* - 

IVl^g — iWg4 — U, ^^410 — ^ -j- 


72 


72 ’ 


Ml = Ml = %{v^A: - yv2), Ml = Ml = hlh'l 


72 ' 


^58 ~ "^85 ~ 0 ) "^59 ~ "^95 ~ ^3 ^ 6 ) ^510 ~ "^105 ~ 0 ) ^67 ~ ^ 7 % ~ 0 ) 

Ml = Ml = h',K*, Ml = Ml = 0, Mlg = Mil = Khl Ml = Ml = ^v^A*, - yv^) . 


72' 


Ml = Ml = hlh, Mlg = Ml, = 0 

Mgg = Mgg = 0, Mg 4 Q = Mggg = k^kj, Mg,Q = M,gg = -^(^1744^ — /iW2) • 

We can diagonalize this hermitian mass squared matrix by the unitary transformation 

= diagiMl.M'l, M|, M|, M|, M|, M|, M|, M|, M?^) , (89) 

Next we write the mass^ matrix in the sups sector the basis (7, Tl, in, Tr, cl, c_r, ul, ur, Ul, Ir)- 
Thus here we denote the sups mass^ matrix in the form {M?)ij = rnfj where 
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,2 

^22 


,2 

^33 


044 


n Vn 

/\/f2 _|_2 

\y[\ 

iL + 

?;2| 

M? + ^ 

2 

1/2 P 

n Wo 

/\/f2 j_^ 

NJ 

bO 

2 V? 

/\/f2 _|_Ij 

2 

I/ 2 P 

Q + 

2 


+ |/i 3 |^ + m|cos2/3 ( - - -siia^Ow 


+ |/i5| + 1 /^ 5 ! + 1 ^ 5 ! + l^sl -COS2/3sin 


+ l^ 5 p + cos2/3sin^ 6*^^, 
o 


+ 1 ^ 3 ! + \K\ + \h'^\ + |/i6| - cos 2/3 ( - - - sin 


vlWs? 

2 


"^55 = ^2L + 


"^66 = ^i-2 + 


+ l/igp + m|cos2/3 Q - ^sin^^^^ , 


2 + l/i'sP +-m|cos2/3sin2 0H/, 


rriyy = M|^ + + l/igp + m| cos 2 ^ (- — - sin^ 9 

2 \ 2 3 


IV 


mlf, = Ml + ^^ 2 ^ + l/i'sT + ^^1 cos 2/3 sin^ 9w, 

1 2 
2 ~ 3 


vlW,? 


^99 = ^II + ^ o + l^6p + ^1 COS 2/3 ( - - - sin^ 0 


‘4L ■ 2 

,„2L,/ |2 

‘*4 " 


vy 


2 j \^2 I I I!, |2 I 2 2 0/3 • 2/1 

= - 2 -^ l^sl +-m^cos2/3sin 6 'h/. 
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2 

2* 

77712 

= 77121 = • 

2 

2* 

77715 

= 77751 = / 

2 


77723 

~ 77132 — ( 

2 


77726 

= 77702 = ( 

2 

2* 

77734 

= 77743 = 

2 


77737 

— 77773 ~ * 

2 

2* / 

77745 

= 77754 = ( 

2 


77747 

= 77774 = ( 

2 

2* 

77756 

= 77705 = 

2 

2* 

77757 

= 77775 = ^ 

2 


77767 

— 77770 ~ * 

2 

2* 

77778 

= 77787 = 

2 

2* 

77719 

= 77794 =; 

2 

2* 

77729 

= 77792 = ■ 


2* 

77210 

= 777192 = 


2* 

77310 

= 777493 = 

2 


77749 

— 777g4 — ( 

2 

2* 

77759 

= 77795 =; 

2 


77769 

~ 777 g0 — ( 

2 

2* 

77779 

= 777g7 = ^ 

2 

2* / 

77789 

= 777g8 = ( 

2 

2* 

77910 

= 777i0g = 


viy2hl V2h5y[* ^ o* v'l / , _ x _2 




+ = mfi = = 0, 


^Tg = = 0> "^?7 = = hshl, mil = ^sl = 0> 

n 2 2* y2 / A* ^ 2 2* '^ 1 ^ 3?/2 , "^ 21 / 3^5 

0, m24 = m42 = - m,m25 = = - ^ , 


2 _ 2 * _ '^1^3?/2 I "^ 21 / 4^5 2 


= 0, m27 = 7717 ; = 


72 


+ 


, "^28 = "^82 = 0, 


72 


^1^51/2 _ V_2]m ^2 _ ^2* _ Q ^2 _ ^2. _ u^'* 

72 72 ’ ~ ~ ~ ~ ’ 

0,ml^ = mil = h^K*, 

n 2 2 * y'iv 2 hl , viy 2 h'^ 

0, m46 = 77104 =-^ + 


r\ 2 2* 

= 0, 77748 = 77784 = 


72 

viy2h'i* 

72 


V2y'ih'' 


72 ’ 

// 

3 


72 ’ 


^(7)274* -/iVl), 

^3^3 ) ''^58 ~ ''^85 ~ 0) 
0) ^68 ~ ''^86 ~ ^5^5 ) 

^{v 2AI- y,Vi), 

= h^h^^miiQ = 777^04 = 0, 

7/^7;i/76 ^27/5^8 


72 


72 ’ 


^5'^8’ 

n 2 _ 2 * _ y 5 V 2 h 6 Viy 2 hl 

0 , 77744 Q — 777494 — ^ + 


72 

^6^3 5 77754Q = 777495 ~ 0 
0)777049 = 777490 ~ ^5^8 
/76 7 '* ,777749 = 777497 = 0) 
0)777849 = 777498 ~ ^5^8) 
^(^2^*4 -/iVi) 


72 ’ 
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We can diagonalize the sneutrino mass square matrix by the unitary transformation 




diag(M„? ,M? , Ml, Ml, Ml, Ml, M? , ) 


(91) 
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